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ABSTRACT 

We present the local H II region metallicity near the site of the recently discovered multiply lensed supernova 
(SN; “SN Refsdal’’) at redshift 1.49. “SN Refsdal” is located at the outer spiral arm (~7 kpc) of the lensed 
host galaxy, which we have previously reported to exhibit a steep negative galactocentric metallicity gradient. 
Based on our updated near-infrared integral field spectroscopic data, the gas-phase metallicity averaged in an 
intrinsic radius of ~ 550 pc surrounding an HII region ~ 200 pc away from the SN site is 12 + log(0/H)ppo4N2 
< 8.67. The metallicity averaged over nine H II regions at similar galactocentric distances (~5-7 kpc) as “SN 
Refsdal” is constrained to be 12 + log(0/H)ppo4N2 < 8.11. Given the fortuitous discovery of “SN Refsdal” 
in an advantageously lensed face-on spiral, this is the first observational constraint on the local metallicity 
environment of an SN site at redshift z > 1. 

Subject headings: galaxies: abundances — supernovae: general — gravitational lensing: strong 


1. INTRODUCTION 

Metallicity at the site of the supernova (SN) explosion pro¬ 
vides important clues to the physical origins of the SN(e). Re¬ 
ciprocally, the evolution and feedback of SNe provide crucial 
recipes for the chem ical enrichment models of galaxies across 
cosm i c time (e.g.. IScannapieco et all 120061 : Kobayas hi et all 
[20071: iFew et al. 112014. 

Generally, SNe can be attributed to either the thermonu¬ 
clear ex plosion of a white dwarf in a binary system (Type 
la; e.g., Hillebrandt & Niemeyer 2000: IParrent et alJ 120141) 
or the c ore collapse (CC) of a massive star (Type II , 
Ib/c; e.g.. iHeger et al J [20031 ICrowtherl 120071: ISmarttfl2009f) . 
Theoretically, at low metallicities ([Fe/H] < —1.1), SNe 
la may be suppressed through the single-degenerate chan¬ 
nel as the white dwarf wind would be too weak for the 
mass to reach the Chandrasekhar limit (Kobayas hi et all 
[19981: iKobavashi & Nomotol l2009h . For CC SNe, the 
number ratio of SNe Ib/c to SNe II is expected to in¬ 
crease with progenitor metallicity due to the dependence of 


mass-loss rate on metallicity (e.g.. 

Mevnet & Maedeil 2003; 

Prantzos & Boissier! 2003| Eldridge 

2007; Limongi & Chieffi 

2009, see also Eldridge et al. (2013)); Yoon et al. (2010) for 


the binary scenario). 


However, the strong theoretical dependence on metallic¬ 
ity has not been confirmed in observations. At low redshift 
(z < 1), observations have shown that the metallicity distribu¬ 
tions of the host galaxies of SNe II and SNe la are statistically 
indistinguishable, whe reas the host galaxies of SN Ib/c favor 
higher met allicity (e.g JPrieto et al.ll2008l : lModiaz et all201 lb 
IStoll et akll2013l) . The class of abnormally luminous SNe Il/Ia 
favor lower-metallicity environments (e.g.. IChildress et all 

1201 ltlKhan et alJl201 ltlKellv & Kirshnerll201 2l) . At highred- 

shift (z > 1), measurements for SN host metallicities are 
still rare and deserve special attention (e.g. .iFrederiksen et all 

120141 iGraur et alJl2014tlRodnev et alJl2014lh 

Note that theoretical models of SNe are based on the metal¬ 
licity of the progenitor star, whereas in observations, direct 
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measurement for the progenitor star metallicity is only pos¬ 
sible for SNe in nearby galaxies. Most studies are based 
on assumptions that the global metallicity of the host galaxy 
is a good proxy for the SNe progenitor metallicity and that 
gas-phase oxygen abundance is a good proxy for the stellar 
metallicity. Additionally, the majority of local spiral galaxies 
show a radial metallicity gradient and their gradients could be 
significantly stee per at high redshift (e.g.. [Jones et alJl2010fc 
I Yuan et akll201 lh . It is therefore quantitatively unclear how 
much systematic uncertainty in the observed metallicity and 
SN-type relation is introduced by the approximation of the lo¬ 
cally and globally averaged metallicity (e.g.. IStanishev et ail 
120121 ) . 

In the era of integral field spectroscopic (IFS) surveys, spa¬ 
tially resolved properties at the sites of SNe can now be mea¬ 
sured for large samples at low redshift, enabling more accu¬ 
rate and unbiased local environment parameters to be probed 
for the SN progenitors dGalbanv et al.ll20L4h . Measuring spa¬ 
tially resolved metallicities at high redshift remains challeng¬ 
ing because of the limited a ngular resolution and signal-to- 
noise (S/N) (S/N; [Yuan et aLir2013i) . With a boost of a few 
tens in flux and area through lensing magnification, gravita- 
tionally lensed targets are the forer unners of this field (e.g., 
Hones et aPl2013l : IWuvts et alll2014l) . 

The recent discovery of a m ultiply lensed supernova at 
z= 1.49 (“SN Refsdal”: [Kelly et~akl 120151) in the cluster field 
of MACS J1 149. 5+2223 provides a rare opportunity to study 
the local metallicity environment of a high-redshift SN for 
the first time. In this Letter, we update the IFS data from 
lYuan et akl ( 2011 ) and place local metallicity constraints near 
the explosion site of “SN Refsdal.” Throughout this paper, 
we adopt a flat cosmology with Hm= 0.3, Ha=0.7, and 
Ho=70 km s -1 Mpc~_^ We use solar oxygen abundance 12 
+ log(0/H) 0 =8.69 ( Asplnnd et al1l2009h . 

2. NEAR-INFRARED IFS DATA 

Our data were collected on the laser guide star adaptive 
optics (LGSAO) aided near-infra red integral field spectro¬ 
graph OSIRIS (iLarkin et alJl2006l) at KECK II in 2010 and 
2011. The data presented in lYuan et alJ (1201 lh were taken 
on 2010 March 3. We obtained a total of 4.75 hr net expo¬ 
sure in the Hn3 band on the lensed spiral galaxy “Spll49” 
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Figure 1. Left: OSIRIS Ho; (red) and H ST F814W (green), F555W (blue) band color image of the host galaxy of “SN Refsdal” at z = 1.49 behind the 
lensing cluster MACS J1149.5+2223. The white box shows the 4.8 by 6.4 arcsec FOV of the OSIRIS Hn3 filter with the Of! I per pixel plate scale. The 
foreground elliptical galaxy is masked out (black square). Locations of the “SN Refsdal” images are highlighted with green circles and are labeled as SI, S2, 
S3, and S4. Two Hll regions close to S2 and S3 are labeled as “HlI_SN3” and “HlI_SN2.” Five Ha contour levels (0.8, 1.4, 2.2, 2.8, 3.5, in units of 10“ 16 
ergs s —1 cm -2 arcsec -2 ) are shown in pink. One of the multiply lensed Hll regions defined as associated with the supernova iSharon & Johnsonll2015D is 
labeled as “35.3.” Cyan circles indicate the nine Hll regions that are at approximately similar galactocentric distances (~5-7 kpc) as “SN Refsdal.” A radius of 
Of !3 (~ 550 pc on the source plane) is used to define the HII region size. Right: averaged spectra for different regions marked on the left. Raw data are shown in 
black; the la noise spectra are in red. Cyan lines show the best-fit Gaussian profiles for lines detected above 3 a (S/N of the Ha and [NII] lines are marked). 


((<* 2000 ) 62000 ) = (11 h 49 m 35.284 s , +22°23'45.86")) that is 
hosting “SN Refsdal.” On April 22,2011, we acquired another 
2 hr of net exposure on this source with the same instrumental 
settings under similar weather conditions. The data were re¬ 
duced using the same pipeline and analysis code as described 
in [Yuan et al.l (1201 ll) . We coadded the flux-calibrated data of 
2010 and 2011 and weighted by the S/N of the total Ha flux. 

To summarize, the data presented in this paper are based 
on a total of 6.75 hr on-source exposure on OSIRIS KECK II. 
The spatial resolution in our observation is ~0'/l, correspond¬ 
ing to a source-plane resolution of ~170 pc after correct¬ 
ing for lensing magnification of fj,fi ux =22 (using the same 
lens model as described in Yuan et al. 2011). The spec¬ 
tral coverage is 1.594— 1.676/rm with a spectral resolution of 
~ 3400. The corresponding rest-frame wavelength range is 
6406—6731 A at the redshift of “Spll49” (z=1.49), aiming 
to capture the Ha and [NII] emission lines. 

Figure Q] shows the two-dimensional (2D) Ha map from 
the OSIRIS datacube for Ha line detections of S/N> 3er. We 
use the coordinate s of the SN sites on the image plane from 
iKellv et aD (120151) and lOguril (120151) . The sites of the lensed 
supernova images are marked with green circles (“SI-4”). 
Two SN locations (“S2” and “S3”) are covered in the OSIRIS 


field of view (FOV). We detect a significant amount (/Ha 
> 0.8 x 10 -16 ergs s -1 cm 2 arcsec -2 or > 0.8 M 0 yr -1 
kpc -1 ) of Ha emission from the Hll regions (“Hll_SN2” 
and “Hll_SN3”) that are close to “S2” and “S3.” Most of 
the bright Ha knots are associated with a strong continuum, 
though not all continuum peaks have bright Ha emissions 

(Figure [TJ. 

Note that “H II_SN2” and “H II_SN3” are not the multi- 
ply lensed H II reg i ons as sociated with the SN as defined in 
ISharon & Johnson! d2015l) . We detect very weak (3a) Ha 
emission from one of the multiply l ensed H II regions ( labele d 
as “35.3” in Figure |T]» defined in ISharon & Johnson! (120151) . 
The spectra extracted at the location of “35.3” using the same 
aperture size as “H II_SN2” is similar to “H II_SN2” (see Fig- 
ure[T] panel (b)). 

3. METALLICTTY 

3.1. Oxygen Abundance (12 + log(OZH)) 

We use the [N Il]/Ha ratio and the N2 metallicity calibra¬ 
tion to calculate oxygen abundance (PP04N2; linear relation 
of lPettini & Pagelll2Q04l) . 

In this work, we extract averaged spectra for the HII regions 
(“H II_SN3” and “H II_SN2” on Figure QJ near the SN site. 
We use the peak of the Ha emission as the extraction center. 
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Figure 2. Velocity and velocity dispersion maps derived from Ha lines from the OSIRS datacube. Two Hll regions that are close the SN sites are indicated 


by cyan circles. There is no abnormal velocity structure nor enhanced velocity 
“H IIJSN2” 

A radius of 0f.'3 (— 550 pc on the source plane) is chosen as 
the size of the H II regions empirically because the S/N of the 
Ho line drops below three outside the radius. We test different 
extraction apertures, shapes, and sizes and find no statistically 
significant differences in the averaged spectra if the extracting 
aperture is centered at the peak of Her clumps. 

The 1 a spectra in Figure Q] are derived from bootstrapping 
(50000 times; weighted by Her flux) the individual spectrum 
of each pixel in the extraction aperture. The Gaussian line 
profile fitting was conducted using a ;y 2 minimization proce¬ 
dure with 1/cr 2 as the weighting. We then use the statistical 
error output from the line fitting procedure to estimate the S/N 
of the flux. For line fluxes of S/N < 3, we derive the flux upper 
limit by using the la spectra and median Gaussian line width 
(5.27 A) of the whole galaxy. The statistical metallicity un¬ 
certainties are calculated by propagating the flux errors of the 
[NII] and Her lines. The metallicity calibration of the PP04N2 
meth od itself has ala dispersion of 0.18 dex (Pettin i & Pagell 
120041 ). Therefore, for individual galaxies that have statisti¬ 
cal metallicity uncertainties of less than 0.18 dex, we use the 
systematic errors of 0.18 dex. 

In the right panel of Figure [Q we show the averaged spec¬ 
tra for the H II regions close to the multiple images of “SN 
Refsdal.” Unfortunately, the [NII] line is below the 3a detec¬ 
tion limit for all of the individual H II regions. The projected 
distances from the H II regions to the SN sites are ~ 200 pc 
and 1.2 kpc for “Hll_SN2 (panel (b))” and “Hll_SN3 (panel 
(a))” respectively. “Hll_SN3” is less likely to be physically 
associated with the SN. 

To further improve the S/N on the weak [NII] line, we ex¬ 
tract a spectrum averaged across nine H II regions that are 
at similar galactocentric distances as “H II_SN2,3” (the pro¬ 
jected distance is —5-7 kpc on the source-plane; panel (c) in 
Figure Q}. No [NII] is detected above an S/N of three in 
the averaged spectra across the H II regions, consistent with 
the result from lYuan et alJ ( 2011 ) that [N II |/Hn ratios are ex¬ 
tremely low at the outskirts of the spiral. For comparison, the 
average spectrum from the innermost annulus of ~ 1.8 kpc 
(source-plane) is plotted at the bottom right panel of Figurejl] 


dispersion indicative of strong shocks/winds at the locations of “H II.SN3” and 

We note that analysis on the source plane gives consistent 
results though the S/N of the spectra is slightly downgraded 
by 1.2 x due to further uncertainties in lensing models. We 
conclude that results based on the emission line ratios are not 
affected by the lensing models. 

The absolute value of oxygen abundance depends on the 
calibration method dKewlev & Ellisorj 12008 ). The advan¬ 
tages of using PP04N2 for absolute oxygen abundance in¬ 
clude the following. 1) The [Nll]/Ha ratio is not affected 
by extinction, and in most cases, is the only line ratio avail¬ 
able for high redshift observations. 2) The absolute value of 
PP04N2 metallicity is actually closer to direct electron tem¬ 
perature method (Te) m etallicity than other co mmonly used 
strong line diagnostics dKewlev & EllisorJ 12008). The “Te” 
method metallicity is considered to be a more “direct” method 
to calculating the absolute metallicity of galaxies, though it 
could be underes tima ting the value at very high metallici- 
ties (e.g jNicholls et aPl2013h . The caveats of using PP04N2 
to diagnose the absolute oxygen abundance include the fol¬ 
lowing. 1) The [N II |/Ho ratio is affected by ionization pa¬ 
rameters and the presence of shocks. At a fixed metallicity, 
a higher ionization parameter produces lower [N Il]/Ha: ra¬ 
tios, and the existence of shock will elevate [Nil 1/Ho ratios 
(e.g..lMorales-Luis et al .120 1 3tlRich et alJ201lUKewlev et al.1 

120131) . 2) The N/O ratio may vary across galaxies, and 

nitrogen fails to be a good proxy for oxygen abundance 
when the primary nitrog en production is dominating (e.g., 
iMorales-Luis et ah/2013 ). The N2 index should not be used 
at very low oxygen abundance because the above caveats 
are most severe at low metallicities (12 + log(0/H)< 8.0; 
IMorales-Luis et all2013l) . 

Keeping the caveats in mind, the upper limit we derived 
for the SN site is still useful as a guide. The Hll region 
(Hll_SN2) that is ~ 200 pc away from the SN has an up¬ 
per limit metallicity of 12 + log(0/H)< 8.67. Furthermore, 
the upper limit of 12 + log(0/H)< 8.11 from the nine H II re¬ 
gion average implies that the oxygen abundance is extremely 
low: more than —0.5 dex lower than the central 1.8 kpc and — 
0.3 dex lower than the average metallicity of the whole galaxy. 
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We can not quantify the effect of shock or variations in ioniza¬ 
tion parameters. However, we do not find enhanced velocity 
dispersion indicative of strong shocks/winds at the locations 
of “H II_SN3” and “H II.SN2” (see Figured}. 



Figure 3. Histogram and cumulative distribution of [Fe/H ] for the SNe la. 
II, Ib/c sample from IPrieto et all 1200 Bl) and SNe II from IStoll et all 12013 !). 
These local samples show that SNe Ib/c occur in higher-metallicity host 
galaxies than SNe II, whereas there is no statistically significant difference 
between SNe la and SNe II. The KS probability tests of two host metal- 
licity samples being drawn from the same distribution are 5% (Tl-Ib/c), 3% 
(Ia-Ib/c), and 56% (H-Ia). The vertical dashed-dotted line highlights the the¬ 
oretical [Fe/H] limit of <-1.1 below which the single-degenerate path SN la 
is forbidden for M \y />=1 Mm IKobavashi etaD 11998 ). 


3.2. Iron Abundance [Fe/H] 

Because it is iron rather than oxygen that provides the dom¬ 
inant opacity for the white dwarf winds in the late stage 
of stellar evolution (IKobavashi et al.|[~1998l) . we infer stellar 
iron abundance ([Fe/H]) using gas-phase oxygen abundance 
(12 + log(0/H)). 

The [O/Fe] ratios have a well-known relation with [Fe/H] in 
the solar neighborhood, where [O/Fe] is constan t at [Fe/H] < 
—1 and linearly decreases at [Fe/H] > — 1 (e.g. JCavrel et al.l 
120041) . The plateau [Fe/H] is determined from the IMF 
weighted yields of CC SNe (SNII and Ibc) and is 0.45, con¬ 
sistent with observations of stars with 3D/NLTE corrections 
(Kobayashi et al. 2006). With this constant, the upper limit 
of [Fe/H] of the near the site of ”SN Refsdal” is —1.06 for 
12 + log(0/H)< 8.11 (nine HII average). 

We use the prescriptions in IStoll et afl (12013l) for the con¬ 
version: [Fe/H] = -11.2 + 1.25(12 + log(O/H)) PP04 N 2 - 
We list the converted [Fe/H] values for SN sites in Table 
1. Here, we adopt the solar metallicity 12 + log(0/H) = 8.69 
(lAsnlund et al.[2009l) . If we use other values of solar metallic¬ 
ity, e.g ., 12 + log(Q/H) = 8.86 from lDelahave & Pinsonneaultl 
(120061) , the [Fe/H] ratio is ~ 0.2 dex lower. This conversion is 
based on fitting to the linear correlation between [Fe/H] and 
[O/H], which is valid above [Fe/H] = — 1. 

Note that oxygen is one of the a elements (O, Mg, Si, 
S, and Ca), and observed [o/Fe] ratios do depend on the 
system; [o/Fe] is higher at the given [Fe/H] in the Galac¬ 
tic bulge and thick disk, while [o/Fe] is lower at the given 
[Fe/H] in Magellanic Clouds and dwarf spheroidal galaxies 


(e.g. IKobavashi et~akll2006l) . Damped Ly a systems, which 
may correspond to the outskirts of spiral galaxies at high 
redshifts, show low [o/Fe] ratios, although this depends on 
the dust depletion correction (e.g.. IWolfe et al.|[2005 ). This 
[O/Fe] variation is not taken into account in the column 5 of 
Table 1. It is unknown if the relation between O and Fe stel¬ 
lar abundances can be used to convert the gas abundance. We 
therefore caution the limitation of this method. 

4. DISCUSSION 

Current observational studies on local galaxies have shown 
some patterns between the types of SN progenitor and the host 
galaxy metallicities. For example, by directly comparing the 
central oxygen abundance of the h osts of SNe Ib/c and SNe 
la with SNe II in SDSS galaxies, iPrieto et aD (120081) found 
strong evidence that SNe Ib/c occur in higher-metallicity host 
galaxies than SNe II, whereas there is no statistically signif¬ 
icant difference between SNe la and SNe II. Similar con¬ 
clusions are also found in other studies that have metallic - 
ities directly me asured at the SN sites (Modiaz et alJ 120081 : 
IStoll et al.l2013l) . Although normal SN la have not been found 
at [Fe/H] < —1; the metallicity of one of the most metal- 
poor SN Ia host galaxies is 12 + log(O/H)=8.01 (with the R23 
method /Childress et al.ll2011l) . 

To illustrate this, in Figure[3]we plot the histogram and cu¬ 
mulative distri buti on of [Fe/H] for the SNe la, II, Ib/c sam¬ 
ple from [Prieto et ahl (120081) . We derive the oxygen abun¬ 
dance using PP04N2 and convert it to [Fe/H] using the same 
methods described in Section 3.2. We also plot the sample 
of lStoll et al. ( 201 3s) for Type II SN progenitor regions from 
a single-source survey that minimized the host galaxy selec¬ 
tion bias. The metallicity distributions of SNe II from the two 
samples agree well, and there is no distinguishable difference 
between Types la and II. Type Ib/c show a preference toward 
higher [Fe/H]. 

The vertical dashed-dotted line in Figure 0 marks the the¬ 
oretical [Fe/H] limit of 1.1, b elow which the single - 
degenerate path SN la is forbidden (Kobav ashi et alJ ll998). 
Observations have shown evidence tha t suc h a boundary for 
SNe la exists. For example, [Rodney et all ( 2014 ) found a 
low rate of SNe la at z >1 which could be explained by 
a population of single-degenerate-dominated Type la at high 
redshift. Note that certain types of SNe may occur only 
for low-metallicity pr ogenitor stars (e.g. jLanger e t al. 2007). 
IChildress et alJ (1201 ll) found super-Chandrasekhar candidates 
in metal-poor environments. 

The inferred [Fe/H] upper limit of < —1 using the 

nine HII region average is close to the boundaries of the 
IKobavashi et all (119981) limit, causing mild tension for “SN 
Refsdal” to be a normal Type la. The [Fe/H] of “SN Refs¬ 
dal” is at the low end tail of the [Fe/H] distribution of local 
SNe. Given that Type Ib/c usually favor high-metallicity en¬ 
vironments, it could imply that “SN Refsdal” is less likely to 
be a Ib/c in a relative sense. However, the SN type of “SN 
Refsdal” is still unknown. Future follow-up observations will 
determine the SN type of ‘SN Refsdal” and thus lead to more 
meaningful discussion on the metallicity constraints. 

5. CONCLUSION 

We use our updated near-infrared IFS data to place con¬ 
straint on the local metallicity near the site of “SN Refsdal.” 
Based on our updated near-infrared IFS data, the gas-phase 
metallicity averaged in an intrinsic radius of ~ 550 pc sur¬ 
rounding an H II region ~ 200 pc away from the SN site is 12 
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Table 1 

Metallicity and Observed Line Fluxes 


Region 

Flux Ha 

Flux [NII] 

12+log(0/H )PP04N 2 

[Fe/H] 

(i) 

(2) 

(3) 

(4) 

(5) 






Hll SN3Q.2 kpc from SN a ) 

H II SN2 (200 pc from SN) 

H II SN2,3, H II xl-7 (5-7kpc galactocentric) 
Central 1.8kpc 

Whole galaxy c 

0.026T0.004 

0.010±0.002 

0.022±0.001 

0.025T0.002 

0.024±0.0008 

<0.004 

<0.004 

<0.0009 

0.011±0.0012 
0.004±0.001 

<8.45 

<8.67 

<8.11 

8.69±0.05 b 

8.43±0.05 b 

<-0.64 
< -0.36 
<-1.06 
-0.34 
-0.66 


Note. — (1) Region where the spectra are averaged, labeled in Figure 1. (2), (3) Averaged observed flux density over respective 
spatial regions; in units of 10 —16 ergs s —1 cm -2 per spaxel ((//1 2 ). lcr upper limit is derived for line detection below an S/N of 5. (4) 
Metallicity based on the PP04N2 method. (5) [Fe/H] inferred from oxygen abundance (see Section 3.2) a Intrinsic distance corrected 
for lensing, uncorrected for orientation. Typical uncertainty of the distance from lens model is ~ 20%. b This is the statistical error. 
Systematic error from the N2-dex method is ~ 0.18 dex. c The average flux density of the whole galaxy is derived using the same 
method as the HII regions, wi th an aperture size tha t includes the whole galaxy. The total SFR of the galaxy after lensing and dust 
correction is ~ 1.2 Mq yr —1 (Livermore et al|2015f) . 


+ log(0/H)ppo4N2 < 8.67. The metallicity averaged over nine 
H II regions at similar galactocentric distances (~5-7 kpc) as 
“SN Refsdal” is constrained to be 12 + log(0/H)ppo4N2 < 
8 . 11 . 

Given the fortuitous discovery of “SN Refsdal” in an ad¬ 
vantageously lensed face-on spiral, this is by far the best ob¬ 
servational constraint we can set on the local metallicity en¬ 
vironment of an SN site at redshift z > 1. The next genera¬ 
tion of 30 m telescopes, such as TMT and GMT, are ideal for 
pursuing the H II region metallicity of the host galaxy of ‘SN 
Refsdal.” Future spectroscopic and imaging follow-up work 
is required to conclude the supernova type of “SN Refsdal.” 
Our metallicity constraints provide an important parameter to 
understand the physical origin of “SN Refsdal.” 
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icantly improved this paper. T.Y. thanks the useful discussion 
with Mike Childress, Fang Yuan, Richard Scalzo, and D.-C. 
Nicholls. L.K. acknowledges the support of an ARC Future 
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Observatory. The authors wish to recognize and acknowledge 
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